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Abstract. 

The new technique, multi-frequency imaging ( MFI) is developed. In VLBI, Multi-Frequency Imaging (MFI) consists of multi- 
frequency synthesis (MFS) and multi-frequency analysis (MFA) of the VLBI data obtained from observations on various fre- 
quencies. A set of linear deconvolution MFI algorithms is described. The algorithms make it possible to obtain high quality 
images interpolated on any given frequency inside any given bandwidth, and to derive reliable estimates of spectral indexes for 
radio sources with continuum spectrum. Thus MFI approach makes it is possible not only to improve the quality and fidelity of 
the images and also essentially to derive the morphology of the observed radio sources. 



1. Introduction 

The new technique, multi-frequency imaging ( MFI), is a 
powerful tool not only for multi-frequency synthesis (MFS) 
(Conway et al. (1990l), but also for multi-frequency analysis 
(MFA) of the VLBI data obtained from observations on vari- 
ous frequencies. This tool allows us to obtain both high qual- 
ity images interpolated on a given reference frequency inside 
of a given bandwidth, as well as reliable estimates of spectral 
indexes for radio sources with continuum spectrum. This ap- 
proach is very important not only for future Space VLBI mis- 
sions like Radioastron, but also for ground-based VLBI arrays 
like the EVLA and VLBA. MFI approach makes it is possible 
not only to improve the quality and fidelity of the images and 
also essentially to derive the morphology of the observed ra- 
dio sources. At the same time, this approach shows that MFI 
will provide the highest angular resolution possible for a given 
Space VLBI mission. 



2. Statement of the problem 

Let us consider a linear model for intensity Ik pq = I(x p ,y q , Vk) 
of the radio source in a point (x p ,y q ) on the observational fre- 
quency Vk'. 



where & = In (1 + B k ) * J3 k , 

and thus the spectral indexes a pq = a(x p ,y q ^ can be ob- 
tained as 

(h) pq = a pq ■ Qo) pq 

Let us consider a target function 
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where, w k nm = w(u„, v m , v k ) > are weights, Vknm, Vfown is a 
measured and a model visibility function respectively, 
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where, Ak is a gain coefficient for k-th antenna, 
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tp pq is a normalized beam, A is a grid step. 
The problem of the optimization can be presented as a so- 
lution of the following system of linear equations: 
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where v is reference frequency corresponding to the intensity for a yector of intensity (?) = ^ > (f,) , . . . , (j^) ) 7 , 
(h)) p , r " -" 



If the intensity Ik pq in the point (x p , y q ) can be approximated 
by power law as 



where the m-th residual map (D m ) pq can be defined as: 
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then we can present it as 

hpq = (h)pq e ?k """ ~ (h) pq • (l + ^pq) 
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where, 



M-\ 



Dkpq = ^ Wknm ■ Vknm • exp fai (u„X p + V m y 9 )} 

is a k-th "dirty" map at the point \x p ,y q \ 

M-l 

Bk,p-i,q-l - ^ WknmX 
n,m=o 

x exp {2tt; \u n (x p - x,) + v m (y 9 - y,)]} 
is a k-th "dirty" beam at the point (x p - Xt,y q - y/). 

3. Solution of the problem 

The solution of the problem can be presented as an iterative 
procedure for a vector (i) : 
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and the residual 
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Here E = [EfA is a positive defined matrix of maximum 
values of weighted "dirty" beams, E,j = (-B/+jJ , i, 7 = 
0, - 1 ; y is a loop gain. The process of the iteration can be 
completed if e (s_I) < e, where s is a given accuracy. Otherwise 
it is necessary to suppose 5 = 5 + 1 and to calculate the next 

.Conditions of the convergence of the algorithm above is 
0<y<2, 1<N <K. 

The developed algorithm is nothing other than the multi- 
frequency linear deconvolution (multi-frequency CLEAN), it- 
self. This procedure described in more detail by Likhachev et 
al. (2003 1. Notice that the developed algorithm allows to syn- 
thesize and analyze of high-quality VLBI images directly from 
the visibility data measured observed on a few frequencies, 
without analyses of the images of the source itself. In case of 
multi-frequency linear deconvolution, it is possible to synthe- 
size an image of a radio source at any intermediate frequency 
inside any given frequency band. Thus, spectral interpolation 
of the image is feasible. This part of the algorithm is carry out 
the synthesis of the image itself. However, the algorithm also 
makes it possible to obtain an estimate of the spectral index 
for a given radio source, i.e., it implements the analysis of 
the image. It is clear that multi-frequency imaging (MFI) will 
provide the highest angular resolution possible for any VLBI 
project due to its improved (u, v)-coverage. 
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4. Implementation of the linear deconvolution 
algorithm 

The algorithm described above was implemented in the 
software, Astro Space Locator (ASL) for Windows 
(http : //platon. asc . rssi . ru/dpd/asl/asl . htmli. 
It was developed by the Laboratory for Mathematical Methods 
of the Astro Space Center (Likhachev et al. (I2003i \ 



